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A  nonadiabatic  sequential-tunneling  model  is  developed  and  applied  to  explore  the  common  origin 
of  the  transient  behavior  of  electrons  in  quantum-well  photodetectors  in  the  presence  of  different 
time-dependent  external  sources,  including  device  temperature,  electric  field,  and  incident  optical 
flux.  For  the  time-dependent  temperature,  a  counterclockwise  hysteresis  loop  in  the  tunneling 
current  as  a  function  of  the  swept  temperature  is  predicted  and  attributed  to  a  blockade  or  an 
enhancement  of  the  sequential  tunneling  of  electrons  between  quantum  wells  by  the 
space-charge-field  effect  when  the  device  temperature  is  swept  up  and  then  down.  When  a 
time-dependent  electric  field  is  applied,  a  dynamical  breakdown  of  the  photodetectors  is  predicted, 
where  the  peak  of  total  current  linearly  increases  with  the  frequency  of  an  ac  electric  field  from  its 
static  value  under  a  dc  field.  This  is  due  to  the  presence  of  an  additional  dielectric  current,  which  is 
proportional  to  the  oscillation  frequency  of  the  ac  electric  field  and  whose  peak  value  becomes 
laiger  than  the  value  of  the  saturated  tunneling-current  peak  in  the  high-frequency  domain.  Under 
the  dynamical-breakdown  condition,  the  quantum-well  photodetectors  behave  just  like  a  unifoim 
dielectric  medium.  In  the  presence  of  a  time-dependent  optical  flux,  an  emission-current  spike  is 
predicted  as  a  result  of  the  dominant  enhancement  of  the  escape  probability  of  electrons  from 
quantum  wells  over  the  loss  of  electron  density  when  an  applied  dc  electric  field  is  small.  The 
experimental  observations  of  the  transient  behavior  of  electrons  in  quantum- well  photodetectors  are 
successfully  reproduced  by  our  numerical  calculations.  [DOI:  10.1063/1.1415760] 


I.  INTRODUCTION 

In  a  recent  article,  Singh  and  Cardimona  reported  a  dis¬ 
covery  of  a  residual  dark  current  in  quantum-well  photode¬ 
tectors  (QWIPs)  when  an  ac  bias  voltage  sweeps  through 
zero.  Later,  Hubbs  et  al ?  observed  a  rolloff  of  the  dynamical 
responsivity  in  QWIPs  when  the  frequency  of  a  chopped 
incident  optical  flux  is  increased  beyond  a  certain  value. 
More  recently,  Choi  et  al?  experimentally  found  a  counter¬ 
clockwise  hysteresis  loop  for  the  tunneling  current  and  a 
clockwise  hysteresis  loop  for  the  emission  current  in  QWIPs 
as  the  device  temperature  was  swept  up  from  1*0  to  300  K 
and  then  back  down. 

A  phenomenological  circuit  model1  was  proposed  that 
successfully  explained  the  observed  zero-bias  residual  dark 
current  in  numerical  simulations.  After  that,  a  quantum- 
mechanical  model4  was  developed  to  explore  the  micro¬ 
scopic  origin  of  the  residual  dark  current.  This  led  to  the 
discovery  of  the  current  instability  and  hysteresis  loop  in 
QWIPs.  For  the  rolloff  of  the  responsivity  in  QWIPs  as  a 
function  of  the  chopper  frequency,  Huang  et  al?  have  put 
forward  a  quantum  theory  that  numerically  reproduced  and 
physically  explained  this  experimental  phenomenon.  A 
scheme  based  on  the  nonadiabatic  model4,5  was  proposed  to 
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eliminate  this  responsivity  rolloff  using  a  small  ac  bias  to 
compensate  the  charge-density  fluctuations  in  quantum  wells 
(QWs)  that  occur  with  a  changing  incoming  optical  flux. 

Here,  we  raise  the  following  question:  Is  there  a  com¬ 
mon  microscopic  origin  for  all  these  experimental  phenom¬ 
ena,  including  the  zero-bias  residual  dark  current,  the  rolloff 
of  the  dynamical  responsivity,  and  the  thermal  hysteresis 
loops  for  both  the  tunneling  and  emission  currents  in 
QWIPs?  We  have  found  the  answer  to  this  question  to  be 
yes,  the  charge-density  fluctuations  in  QWs  induced  by  a 
nonadiabatic  sequential-tunneling  or  photoemission  process. 
It  is  well  known  that  resonant  electron  tunneling  in  QWIPs 
can  occur  only  when  the  barrier  between  adjacent  wells  is 
thin.  If  the  barrier  is  very  thick,  the  phase  of  the  wave  func¬ 
tion  will  be  completely  lost  as  an  electron  tunnels  from  one 
well  to  the  next.  As  a  result,  only  sequential  tunneling  of 
electrons  exists  for  thick  barriers.  If  the  QWIPs  are  only 
subject  to  constant  external  sources,  such  as  device  tempera¬ 
ture,  bias  voltage,  and  incident  optical  flux,  elections  in  the 
QWs  remain  in  a  steady  state  with  a  fixed  electron  density. 
However,  when  the  external  sources  become  time  dependent, 
electrons  in  QWs  can  no  longer  remain  in  this  steady  state. 
Instead,  electrons  will  migrate  through  a  series  of  intermedi¬ 
ate  transient  states.  The  chemical  potential  of  electrons  in 
any  one  of  these  states  suffers  a  fluctuation  from  the  time- 
dependent  external  sources,  and  electron  sequential  tunneling 
between  adjacent  wells  changes  from  an  adiabatic  process  to 
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a  nonadiabatic  one.  As  a  result,  the  charge  density  of  elec¬ 
trons  in  the  QWs  starts  to  fluctuate  around  the  value  of  the 
doping  density,  which  either  blocks  or  promotes  the  tunnel¬ 
ing  of  electrons  via  a  nonadiabatic  space-charge-field  effect 
when  the  QWs  are  charged  or  discharged,  respectively. 

In  this  article,  we  propose  a  nonadiabatic  sequential- 
tunneling  model  containing  a  chemical-potential  fluctuation 
for  QWIPs  in  the  presence  a  time-dependent  external  source. 
On  the  basis  of  this  model,  we  derive  a  dynamical  equation 
with  respect  to  a  nonadiabatic  electric  field  describing  the 
charge-density  fluctuations  in  QWs.  Using  the  derived  dy¬ 
namical  equation  subject  to  a  time-dependent  device  tem¬ 
perature,  we  successfully  reproduce  and  explain  the  observed 
hysteresis  loop3  for  the  tunneling  current  as  a  function  of 
device  temperature  in  QWIPs.  A  completely  different  physi¬ 
cal  mechanism  behind  this  thermal  hysteresis  loop  will  be 
provided  compared  to  the  interpretation  given  by  Choi  et  al? 
We  also  report  on  the  dynamical  breakdown  of  QWIPs  ob¬ 
served  in  our  experiment  and  we  elucidate  the  physics  of  this 
phenomenon.  In  addition,  an  explanation  of  the  previously 
observed  emission-current  spike  by  Hubbs  et  al 2  when  an 
optical  shutter  is  opened  will  be  provided. 

The  organization  of  this  article  is  as  follows.  In  Sec.  II, 
we  present  our  nonadiabatic  sequential-tunneling  model  and 
theory  to  study  the  transient  behavior  of  QWIPs  in  the  pres¬ 
ence  of  a  time-dependent  device  temperature,  electric  field, 
and  incident  optical  flux.  The  thermal  hysteresis  loop,  dy¬ 
namical  breakdown,  and  emission-current  spike  observed  ex¬ 
perimentally  in  QWIPs  are  successfully  reproduced  and  the 
common  origin  of  all  these  transient  behaviors  is  explored. 
Numerical  results  and  discussions  are  given  in  Sec.  HI,  for 
three  different  cases,  including  an  ac  electric  field,  a  chopped 
optical  flux,  and  a  time-dependent  device  temperature.  Some 
previous  and  current  experimental  data  from  other  groups 
and  our  experiments  are  also  displayed  in  this  section  for  a 
qualitative  comparison.  The  article  is  concluded  in  Sec.  IV 
with  some  remarks. 


II.  MODEL  AND  THEORY 

In  this  section,  we  start  by  considering  a  typical  QWIP 
sample  which  consists  of  GaAs  layers  for  the  QWs  and 
A^Ga^As  layers  for  the  barriers  separating  adjacent  wells. 
The  thickness  (width)  of  the  barriers  (wells)  is  LB(LW),  and 
the  depth  of  the  wells  is  U0 .  The  effective  mass  of  electrons 
is  m *.  An  electric  field  is  applied  to  the  sample  in  the  z 
(growth)  direction.  The  donors  are  uniformly  distributed 
within  each  QW  and  assumed  to  be  ionized  to  give  rise  to  an 
electron  density  n2o  in  a  QW. 

The  sequential-tunneling  electrons  in  a  QWIP  will  re¬ 
main  in  a  steady  state  when  a  constant  electric  field  or  a 
constant  optical  flux  has  been  applied  to  the  system  at  any 
fixed  temperature.  In  this  case,  a  steady-state  current  will 
flow  between  the  emitter  and  collector  of  the  system,  and  the 
charge  density  in  each  QW  will  be  kept  at  a  constant  n2D  (the 
applied  electric  field  is  assumed  to  be  uniform  inside  the 
QWIP).  However,  a  small  current  surge  which  is  superposed 
on  this  steady-state  current  will  be  created4,5  in  the  system 
whenever  the  chemical  potential  /jlc  of  electrons  in  the  sys¬ 


tem  undergoes  a  fluctuation  A  /ulc  due  to  an  external  time- 
dependent  source  such  as  device  temperature  Te(t)  or  elec¬ 
tric  field  £b(t ).  The  current  surge  is  so  small  that  its  effect 
cannot  be  seen  in  a  resonant-tunneling  process  within  a  su¬ 
perlattice.  However,  in  a  sequential-tunneling  process  the 
current  surge  can  induce  up  to  20%  charge-density  fluctua¬ 
tions  in  the  QWs,5  which  leads  to  a  measurable  nonadiabatic 
tunneling  current.  By  using  Levine’s  sequential-tunneling 
model6  for  electrons,  the  tunneling-current  surge  can  be 
quantitatively  described  by4 

I  m*  \  f  +00 


X  /« 


E+Ex- 


kBTe 


“/o 


E+Ex-  fic 


(1) 


where  fo(x)  =  (l+ex)~{  is  the  Fermi-Dirac  distribution 
function,  S  is  the  area  for  the  cross  section  of  the  sample, 
vd(£b)  ls  the  drift  velocity  of  electrons4  under  an  electric 
field  £b,  Ex  is  the  ground-state  energy  of  electrons  in  each 
QW,  and  7{E,£b)  is  the  instantaneous  transmission  coeffi¬ 
cient  of  electrons  through  a  biased  barrier  between  adjacent 
wells.7  During  the  sequential-tunneling  process,  an  electron 
sees  an  instantaneous  value  of  the  electric  field  due  to  an 
extremely  fast  tunneling  process  (in  the  range  of  1  ps)  com¬ 
pared  to  the  time  (in  the  range  of  1  ms-1  s)  required  for  the 
appreciable  change  in  £b(t).  In  order  to  see  an  explicit  time 
dependence  of  7 \E,£b),  a  microwave  field  with  the  oscilla¬ 
tion  frequency  exceeding  1-10  MHz  must  be  employed, 
where  a  time-dependent  Schrodinger  equation  must  be 
solved  to  include  the  time  dependence  in  the  transmission 
coefficient  through  a  dynamical  barrier.8 

Equation  (1)  of  our  model  is  based  on  the  following 
facts.  First,  there  exists  a  fluctuation  A/jlc  in  the  chemical 
potential  of  the  electron  gas  in  each  QW  whenever  either  the 
bias  field  or  the  temperature  varies  with  time.  Second,  A  julc 
will  further  introduce  a  conduction-current  surge  Ist .  As  a 
result,  it  causes  a  charge-density  fluctuation  in  each  quantum 
well,  which  in  turn  induces  a  time-dependent  space-charge 
field.  Finally,  the  generation  of  the  dynamical  space-charge 
field  should  be  determined  by  the  current-charge  conserva¬ 
tion  law.  This  model  has  been  successfully  applied  to  explain 
the  zero-bias  offset  of  the  tunneling  current4  and  the  fre¬ 
quency  rolloff  of  the  responsivity5  observed  experimentally. 
For  fixed  electron  density  in  each  QW,  fic  is  a  function  of  the 
electron  temperature.  If  the  temperature  suffers  from  a  fluc¬ 
tuation,  it  induces  A^ac .  Similarly,  if  a  bias  field  suffers  from 
a  fluctuation,  some  electrons  will  be  accelerated  while  others 
will  be  decelerated.  The  imbalance  of  the  electron  energy 
disturbs  the  electron  distribution  and  causes  A/ulc  .  If  the  pe¬ 
riod  of  a  time-dependent  electric  field  is  much  larger  than  the 
tunneling  time  but  smaller  than  the  charging/discharging 
time  of  electrons  in  QWs,  the  electron  motion  is  nonadia¬ 
batic  in  nature.  The  origin  of  the  nonadiabatic  motion  of 
electrons  is  lst  represented  by  Eq.  (1). 

We  would  further  like  to  point  out  that  in  the  Levine’s 
description6  of  the  quantum  sequential-tunneling  current  of 
electrons  the  two-dimensional  tunneling  current  is  effectively 
treated  as  a  three-dimensional  tunneling  density  (3DTD)4 
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moving  with  a  drift  velocity  under  a  bias  voltage  through  a 
thick  barrier  layer  between  two  QWs  in  QWIPs,  as  can  be 
seen  from  Eq.  (1).  The  electrons  contributing  to  3DTD  have 
energies  expanded  from  the  edge  of  the  ground  subband  up 
to  infinity,  which  includes  the  contributions  of  both  thermi¬ 
onic  and  tunneling  currents,  although  the  density-of-states  in 
Levine’s  formula6  can  be  replaced  by  a  3D  one  when  the 
electron  energy  is  above  the  barrier.  The  effect  of  the  energy 
barrier  to  electrons  is  included  as  a  tunneling  transmission 
coefficient  in  3DTD,  while  the  spatial  scattering  of  electrons 
by  imperfections  in  an  alloy  barrier  layer  is  included  in  the 
drift  velocity.  Therefore,  the  drift  velocity  is  appropriate  for 
not  only  the  electrons  above  the  barriers  but  also  the  electron 
within  the  barriers. 

Photoexcited  emission  current  exists  in  QWIPs  as  a  re¬ 
sponse  to  an  incident  optical  flux.  The  generation  of  an  emis¬ 
sion  current  is  probably  a  process  involving  photon-assisted 
electron  tunneling  as  indicated  by  Levine.6  It  consists  of  two 
successive  processes.  The  first  step  is  the  optical  absorption 
by  electrons,  which  lifts  the  electrons  from  a  lower  ground 
state  to  an  upper  excited  state  in  the  QWs.  The  second  step 
involves  the  tunneling  of  photoexcited  electrons  out  of  the 
QWs  through  a  biased  barrier  to  a  continuum  state.  This 
mechanism  for  producing  an  emission  current,  which  is  due 
to  nonadiabatic  photoexcitation,  is  obviously  quite  different 
from  that  for  producing  the  tunneling  currents  described  pre¬ 
viously.  Emission  current  is  excited-state  transport,  while 
tunneling  current  is  transport  in  the  ground  state.  The  de¬ 
tailed  form  of  Eft)  will  be  given  in  Sec.  II C. 

The  existence  of  either  Ist(t)  or  Ise(t)  will  cause  a  devia¬ 
tion  of  the  tunneling  current  in  the  system  away  from  its 
steady-state  value  and,  consequently,  fluctuations  of  the 
charge  density  of  electrons  in  each  QW.4,5  Here,  the  electrons 
no  longer  remain  in  the  initial  steady  state  with  a  time- 
independent  tunneling  current  but  instead  move  through  a 
series  of  intermediate  transient  states  with  charge-density 
fluctuations  and  a  dynamical  tunneling  current.  However,  the 
induced  charge-density  fluctuations  in  the  QWs  will  be  com¬ 
pensated  by  the  concomitant  change  of  the  time-dependent 
tunneling  current  compared  to  its  initial  steady-state  value 
and  the  system  will  eventually  reach  a  final  steady  state  with 
a  different  static  tunneling  current  after  a  characteristic  time 
(in  the  range  of  1  s)  similar  to  the  charging/discharging  time 
of  a  capacitive  system  in  classical  electrodynamics.1 

A.  Time-dependent  device  temperature  Te(t) 

When  the  electron  temperature  Te{t)  is  varied  with  time 
t,  the  chemical  potential  jic(Te)  (for  simplicity,  the  depen¬ 
dence  of  ilc  on  the  fixed  electric  field  is  not  shown)  of  elec¬ 
trons  in  each  QW  is  shaken  up  by  a  fluctuation  A  fic 
~Tr[dpcldTe'\[dTe(t)ldt'],  where  Tr~0.1-1  ps  is  a  param¬ 
eter  reflecting  the  energy-relaxation  time  of  electrons  back  to 
an  equilibrium  state  through  intrasubband  scattering.9  If  a 
parabolic  energy  dispersion  is  adopted  for  the  electron  mo¬ 
tion  within  the  quantum-well  plane  (perpendicular  to  the 
growth  direction),  we  find  the  chemical  potential  of  electrons 
for  fixed  n2D  to  be 


Pc(Te)-Ei+kBTe]n 


Trh2n2D\ 

m*kBTe) 


(2) 


From  Eq.  (2)  we  find  dfic/dTe< 0.  By  using  Eq.  (1),  the 
dynamics  of  the  charge-density  fluctuation,  described  by 
£,yJ  t)  in  each  QW  is  found  to  be  determined  by4,5 
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where  the  dynamical  quantum-well  capacitance  is 


(3) 
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In  Eq.  (3),  £aft)  is  a  nonadiabatic  electric  field  which  is 
associated  with  the  charge-density  fluctuation  in  each  QW,4 
and  /,[£),  ,7’(.(/)]  can  be  calculated  using  Levine’s 
sequential-tunneling  model.4-6 

The  capture  process  of  carriers  flowing  above  the 
barriers10  has  not  been  neglected  in  Eq.  (4)  although  it  has 
not  been  explicitly  shown.  For  the  tunneling  process  [includ¬ 
ing  either  a  time-dependent  temperature  Te(t)  here  or  a  time- 
dependent  electric  field  £b(t)  in  Sec.  II B],  the  difference 
between  the  capture  and  tunnel-emission  currents  is 
E[£b(t)  +  £m({),Tft)]-It[£b(t),Tft)].  Here,  we  would 
like  to  emphasize  that  the  difference  between  the  capture  and 
emission  currents  originates  from  the  nonadiabatic  current 
Surge  in  our  model.  For  the  tunneling  process,  the  current 
surge  results  from  A fic  [<* dTft)/dt  here  or  d£b(t)ldt  in 
Sec.  HB]  in  each  QW.  The  nonadiabatic  current  surge  sur¬ 
vives  even  under  a  uniform  time-dependent  electric  field  and 
it  reflects  the  quantum  property  of  electrons  in  the  QWs.  This 
quantum  effect  is  absent  in  the  previous  theories  under  a 
nonuniform  electric  field,11’12  where  the  change  of  charge 
density  with  time  in  QWs  results  from  the  redistribution  of  a 
nonuniform  bias  field  at  each  time. 

The  simplest  example  of  a  time-dependent  device  tem¬ 
perature  to  help  visualize  the  physics  involved  in  Eq.  (3)  is 
the  step  function:  Tft)  =  T1  +  (T2-T1)6(t-t1 )  which  steps 
at  t=tx.  From  Eq.  (3)  we  know  that  electrons  with  chemical 
potential  /xfTf  are  knocked  away  by  Eft)  at  t=tx  from 
their  initial  steady  state  due  to  dT ft  x)ldt  + 0.  As  a  result,  the 
charge  density  in  each  QW  fluctuates  as  t>tx  due  to  the 
tunneling-current  surge.  If  dTftx)/dt>Q(T2>Tx),  we  find 
£naM<0  for  t>tx ,  which  means  electrons  will  be  sucked  in 
and  leads  to  the  QWs  being  charged.  During  this  time,  the 
electrons  are  no  longer  in  the  initial  steady  state  with  a  time- 
independent  tunneling  current  but  instead  have  moved 
through  a  series  of  intermediate  transient  states.  The  buildup 
of  electrons  in  the  QWs,  however,  decreases  the  tunneling 
current  /,[4+5na(r),r2] </,[£*  ,r2]  with  time  for  t>tx . 
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This  will  bring  in  less  electrons  to  each  QW  and  eventually 
stabilize  the  electrons  to  a  final  steady  state  with  a  higher 
tunneling  current  It[£b  ,T2]  and  different  chemical  potential 
!jlc(T2)  but  the  same  electron  density  «2D.  The  facts  that 
a^+^na(0,T2]</r[^,r2]  with  TX<T2  and  It[£b 
+  5na(0,r2]>/f[4,r2]  with  Tx>T2  during  a  series  of  in¬ 
termediate  transient  states  give  rise  to  a  counterclockwise 
thermal  hysteresis  loop 3  in  It[£b+£ na(0>Te(t)]  versus  Te(t) 
when  the  electron  temperature  Te(t)  is  swept  up  and  then 
down  with  time.  In  this  case,  the  tunneling  current  at  a  cer¬ 
tain  temperature  Te(t)  depends  on  the  history  of  the  variation 
of  Te(t),  namely,  dTe(t)/dt. 

B.  Time-dependent  electric  field  £b(t) 

If  the  applied  electric  field  £b{t)  is  time  dependent,  the 
chemical  potential  p c(£b)  (for  simplicity,  the  dependence  of 
pc  on  the  fixed  electron  temperature  is  not  shown)  of  elec¬ 
trons  in  each  QW  again  fluctuates  with  A fic=  rt[dfjbcl 
d£b][d£b(t)/dt ],  where  rr~0.01-l  ps  is  a  parameter  asso¬ 
ciated  with  the  electron  sequential-tunneling  time,  which  can 
be  evaluated  using  the  ratio  of  electron  charge  to  the  tunnel¬ 
ing  current.13  For  a  parabolic-dispersion  model,  we  find  from 
Eq.  (2)  that  dpc/d£b  =  dE1(£b)ld£b  for  fixed  Te  and  n2D , 
which  relates  to  the  Stark  shift14  of  the  ground-state  energy 
of  electrons  in  the  QWs.  With  the  help  of  Eq.  (1),  the  dy¬ 
namics  of  the  charge-density  fluctuation  in  the  QWs  can  be 
obtained  from4,5 
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where  is^O)  and  pc{ 0)  are  the  ground-state  energy  and 
chemical  potential  of  electrons  in  QWs  at  £b  =  0  and 
can  be  calculated  using  Levine’s  sequential- 
tunneling  model 4-6  If  we  formally  introduce  an  effective 
quantum-well  capacitance  according  to 


y  </[£*(')  ]t/ 

dE  {{£b) 

f  m*eS 

eff 

lblw 

d£b 

[  7 rh2kBTe 

i; 


X  dET[E,£b(t)]f0 


E+Em-^i  0) 


kBTe 


l-/o 


E+EiiOyjt^jQ) 

kBTe 


(6) 


and  notice  that  r,T[£:,£:fr(f)]~Ls/i)ii[ft(0]  and 
dEy{£b)!d£b~eLw,  we  get  Ceff~CQW[0]  as  can  be  ob¬ 
tained  from  Eq.  (4)  by  setting  5na(f)=0.  Under  this  condi¬ 
tion,  Eq.  (5)  will  reduce  to  the  previous  dynamical  equation.4 

We  can  easily  understand  the  physics  implied  by  Eq.  (5) 
again  using  a  step  function  £b(t)  =  £l  +  (£2  —  £l)6(t-tl) 
which  steps  at  time  t—t\.  For  t<t1 ,  the  electrons  remain  in 
an  initial  steady  state  with  a  relative  chemical  potential 


• 

fJ'c(^i)~E{(Ei)  =  fjic(0)  —  Ei(0)  (independent  of  £b(t)  al¬ 
though  the  electron  wave  functions  depend  on  it4)  in  each 
QW.  Because  d£b(t{)/dt^  0,  the  electrons  in  the  QWs  are 
pushed  away  from  their  initial  steady  state.  This  leads  to 
charge  density  fluctuations  in  each  QW  for  t>tx.  If 
d£b(tx)/dt> 0(£2>£x)  is  assumed,  we  find  Em(t)>0  for  t 
>ty ,  indicating  a  discharging  status  for  the  QWs.  However, 
this  will  increase  the  tunneling  current  /f[£2+£na(0>Tj 
>/r[£2  with  time  for  t>tx ,  which  pours  more  electrons 
into  the  discharged  QWs  to  compensate  for  the  charge  loss 
and  helps  the  electrons  establish  a  final  steady  state  with  a 
higher  tunneling  current  It[£2>Te]  but  the  same  relative 
chemical  potential  fic{£2)-Ex{£2)- ilc{0)-Ex{{))  and 
electron  density  n2D .  From  Eq.  (5)  we  know  that  the  tunnel¬ 
ing  current  It[Eb(t)  +  £na(t),Te]  in  the  system  depends  on 
the  nonadiabatic  electric  field  £m(t)  and  the  oscillation  fre¬ 
quency  f 1b  of  £b{t)  through  d£b(t)ldt.  Our  numerical  stud¬ 
ies  have  discovered  that  the  peak  of  It[£b(t)  +  £Jit)  »rj  as  a 
function  of  fl^  saturates  to  a  maximum  value  ~/,[(l 
+  Qfr^Qw[0])^0’Te]  when  ilb  is  very  large,  where  £0  is 
the  amplitude  of  a  sinusoidal  £b(t).  Within  the  saturation 
regime,  the  peaks  of  the  tunneling  current  and  £b(t)  are  in- 
phase  with  each  other  completely,  and  the  maximum  value  of 
the  tunneling  current  changes  with  Te  and  £0.  On  the  other 
hand,  we  know  from  classical  electrodynamics15  that  there 
exists  a  dielectric  current  €0€rSd£b(t)/dt  flowing  through  a 
QWIP  if  the  well  and  barrier  layers  between  the  emitter  and 
collector  are  viewed  as  a  uniform  dielectric  medium  (in  ad¬ 
dition  to  the  QWs  doped  with  electrons)  with  er  being  their 
average  dielectric  constant.  This  dielectric  current  which 
does  not  contribute  to  change  density  in  QWs  increases  lin¬ 
early  with  and  becomes  negligible  when  Elb  is  small.  It 
has  a  rr/2  phase  shift  with  respect  to  £b(t).  Once  the  peak  of 
the  dielectric  current  becomes  larger  than  the  maximum 
value  of  the  tunneling  current  as  is  large  enough,  QWIPs 
suffer  from  a  dynamical  breakdown ,  i.e.,  the  QWBP  will  just 
behave  like  a  uniform  dielectric  medium.  In  this  case,  the 
current  flowing  inside  a  QWIP  is  dominated  by  the  dielectric 
current  instead  of  the  conduction  current. 


C.  Time-dependent  optical  flux  <J>ph(f) 

In  the  presence  of  a  time-independent  electric  field  £b 
and  incident  optical  flux,  the  electrons  in  both  tunneling  and 
emission  channels  of  a  QWIP  remain  in  their  steady  state.  In 
this  case,  the  matching  of  the  emission  and  capture  currents 
keeps  the  electron  density  in  each  well  at  a  constant  n2 D.6 
However,  if  we  apply  a  time-dependent  optical  flux  ®ph(0  to 
the  system,  a  mismatch  between  the  active  emission  and  pas¬ 
sive  capture  currents  occurs5  due  to  a  current  surge,  leading 
to  charge-density  fluctuations.  As  mentioned  at  the  beginning 
of  this  section,  the  mechanism  for  producing  the  emission 
current  is  quite  different  from  that  for  producing  the  tunnel¬ 
ing  current.  Instead  of  using  Eq.  (1),  the  emission-current 
surge  due  to  the  nonadiabatic  photoexcitation  process  is 
found  to  be5 

\  *phW  /  dt 


(7) 
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where  /J&.fcpi/r)]  is  the  steady-state  emission  current, 
which  can  be  calculated  by  using  Levine’s  photoemission 
model,6  and  t,~  0.1-1  ps  is  related  to  the  lifetime  of  elec¬ 
trons  in  the  upper  excited  state.6  lse(t)  will  induce  both  a 
charge-density  fluctuation  in  each  QW  and  a  deviation  of  the 
emission  and  tunneling  currents  from  their  steady-state  val¬ 
ues  Ie[£b  ,d>Ph(t)]  arid  /,[£*,  rj.  By  using  Eq.  (7),  the  dy¬ 
namical  equation  describing  the  charge-density  fluctuation  in 
each  QW  with  dual  conducting  channels  is  found  to  be5 


CQwt^naW] 


d£Jf)_ 

dt 


•  ,rr„  _  ,,/^phW 

Lb* ph(t))  e[£b,^M]  dt 


x{a^,4»ph(0]-/r[^,re] 
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Here,  the  nonadiabatic  emission  current  with  £nd(t)=t0  is 
calculated  by6 
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and  It[£b  ,TJ  is  still  given  by  Levine’s  sequential-tunneling 
model,4”6  where  £i(0)  and  /uc(0)  (for  simplicity,  the  depen¬ 
dence  of  pc  on  the  fixed  electron  temperature  is  not  shown) 
still  represent  the  ground-state  energy  and  chemical  potential 
of  electrons  in  the  QWs  at  £b=0 ,  ftwph  is  the  incident  photon 
energy,  P e[£b]  and  crop[cuph,^]  are  the  escape  probability 
and  optical  cross  section  of  electrons.5 

In  Eq.  (9),  the  capture  process  of  carriers  flowing 
above  the  barriers  has  been  implicitly  included.  In 
fact,  for  the  photoemission  process  with  a  time-dependent 
optical  flux  <J>ph(0,  the  difference  between  the  capture  and 
emission  currents  becomes  /,[£*+ 

+ 1  el ?b  * ' ®  Ph( 0  ]  -  I JL  ?b + £na(  0  *  ®  ph(  0  ]  •  Moreover,  the 
capture  probability  pc  has  been  explicitly  included  in  the 
previous  calculation  for  the  reponsivity  7£sp.5  Here,  we 
would  like  to  indicate  that  the  difference  between  the  capture 
and  emission  currents  in  our  model  comes  from  the  nonadia¬ 
batic  excitation  of  electrons  [<xrf3>ph(r)Afr]  from  the  ground 
state  to  the  first  excited  state  in  the  QW. 

If  we  choose  a  step  function  3>Ph(0  =  ^>i  +  (4>2 
stepping  at  time  t=tu  we  know  from  Eq.  (8) 
that  the  electrons  will  migrate  away  from  the  initial  steady 
state  because  of  d^ph(tx)/dti=0.  Consequently,  the  electron 
charge  density  fluctuates  in  each  QW  as  t>tl.  If 
d$>ph(h)/dt>0(<t>2>&i)  is  taken  as  an  example,  we  find 
0  for  t>tx.  This  implies  that  the  QWs  are  being 
discharged  during  this  time.  The  increasing  tunneling  current 

with  time  for  t>tx  adds  more 
electrons  to  each  QW  to  compensate  for  the  charge  loss. 
Eventually,  the  electrons  in  the  QWs  transit  to  a  final  steady 
state  with  a  higher  emission  current  Ie[£b  ,<J>2]  but  the  same 


TABLE  I.  GaAs/AI0  3Gao  7As  QWIP  sample  1  used  in  our  numerical  calcu¬ 
lation  with  saturation  velocity  vs ,  saturation  field  Ss ,  cross-sectional  area  <S, 
electron  density  h2D ,  and  relative  dielectric  constant  er . 


vs  (106  cm/s) 

£s  (kV/cra)  <S(1CT4  cm-2)  n2D(10ncnTJ)  sr 

2 

2  2.25  4.1  11.2 

chemical  potential  and  electron  density  due  to  the  establish¬ 
ment  of  a  match  between  the  emission  and  capture  currents 
in  the  steady  state.  The  nonadiabatic  emission  current  Ie[£b 
+  ^ia(0»®2]  can  be  either  larger  or  smaller  than  the  steady- 
state  value  of  the  emission  current  Ie[£b  ,<I>2],  depending  on 
the  competition  between  the  loss  of  electron  density  in  the 
QWs,  seen  as  the  term  containing  /*c(0)  -  £L5£na(Q 
~EX( 0)  in  Eq.  (9),  and  the  enhancement  of  the  escape  prob¬ 
ability  /%[£&+ £na(0]  due  to  the  additional  lowering  of  the 
barrier  by  the  nonadiabatic  electric  field.  We  note  that  the 
increase  of  the  nonadiabatic  electric  field  £na(t)  due  to 
d^Ph(ti)/dt>0  makes  the  enhancement  of  the  escape  prob¬ 
ability  dominant  through  its  exponential  dependence  on  it6 
under  certain  conditions.  However,  the  loss  of  charge  density 
in  each  QW,  which  tends  to  reduce  the  final-state  emission 
current  /*[££  + £na(/),<I>2],  will  eventually  dominate  at  a 
later  time  due  to  0^Pe^l.  This  leads  to  an  emission-current 
spike  around  the  stepping  time  tx. 

III.  NUMERICAL  RESULTS  AND  DISCUSSION 

In  our  numerical  calculation,  we  have  chosen  a  QWIP 
sample  1  which  contains  GaAs  layers  as  QWs  and 
Alo.3Gao.7As  layers  as  banders  between  adjacent  wells.  The 
well  depth  is  £/0  =  224  meV,  the  well  width  is  L^=50  A,  the 
barrier  thickness  is  L5=339A,  and  the  effective  mass  of 
electrons  is  ra*  =  0.0665m0,  with  m0  being  the  ffee-electron 
mass.  The  photon  energy  ft<oph  is  resonantly  set  to  be  the 
energy-level  separation  between  the  ground  and  first-excited 
states  in  the  QWs  with  a  homogeneous  level  broadening  of  1 
meV.  The  other  parameters  of  sample  1  are  listed  in  Table  I. 
The  parameter  of  another  sample  used  in  our  measurement 
with  time-dependent  electric  field  is  listed  in  Table  II.  The 
parameter  of  the  sample  used  in  the  experiment  by  Hubbs 
et  al.  with  a  time-dependent  optical  flux  can  be  found  in  Ref. 
2.  In  this  section,  we  will  present  both  our  numerical  and 
experimental  results  for  the  QWIPs  in  the  presence  of  a  time- 
dependent  electric  field.  After  that,  our  calculated  transient 
emission  current  will  be  displayed  together  with  the  experi¬ 
mental  data  from  Ref.  2.  Finally,  we  will  show  our  calculated 
transient  tunneling  currents  with  a  time-dependent  device 
temperature. 


TABLE  H.  GaAs/ Al^Ga^As  QWEP  sample  2  used  in  our  experiments  in 
the  presence  of  a  time-dependent  electric  field  with  well  width  Lw,  barrier 
thickness  LB  ,  cross-sectional  area  <S,  and  electron  density  n2 D . 


Lw  (A) 

L„(  A) 

5(10~4cm-2) 

n2D  (IQ”  cm  2) 

50 

300 

9 

2.5 
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FIG.  1.  Calculated  £na(f)  [left  axis  in  (a)]  and  CQW[£ra(f)]/CQW[0]  [left 
axis  in  (b)]  from  sample  1  at  £0  =  1  kV/cm  and  7^=40  K  as  a  function  of 
Clbt!2ir  for  different  frequencies  £lb/2ir=  0.0025  (thin  solid  curves),  0.025 
(thin  dashed  curves),  and  0.25  Hz  (thin  dash-dotted  curves).  £b(t)  [right  axis 
in  (a)]  is  also  displayed  by  the  thick  solid  curve  for  the  comparison.  Here, 
the  curves  are  only  shown  for  one  time  period  for  clarity. 


A.  Dynamical  breakdown 

For  the  time-dependent  electric  field,  we  choose  a  sinu¬ 
soidal  form 

£b(t)  =  £0sin((lbt),  (10) 

where  £q  and  Clb  are  the  amplitude  and  oscillation  frequency 
of  £b(t).  In  our  numerical  calculation,  the  tunneling  time  rt 
in  Eq.  (5)  is  evaluated  using  e!It[£b(t),Te ],  which  depends 
on  the  time  through  £b(t).  In  the  following,  we  use  sample  1 
in  our  numerical  calculation  and  sample  2  in  our  experimen¬ 
tal  measurement.  In  our  experiment,  an  ac  bias  voltage  is 
applied  to  sample  2  with  its  amplitude  V0  and  oscillation 
frequency  £lb . 

The  solution  of  Eq.  (5)  is  presented  in  Fig.  1,  where  the 
nonadiabatic  electric  field  £na(f)  [left  axis  in  (a)]  and  the 
ratio  of  quantum-well  capacitance  CQW[fna(t)]/CQW[0]  [left 
axis  in  (b)]  are  shown  as  functions  of  £lbt/27r  for  different 
frequencies  ft ^/27T=  0.0025  (thin  solid  curves),  0.025  (thin 
dashed  curves),  and  0.25  Hz  (thin  dash-dotted  curves).  Here, 
we  have  chosen  £0=1  kV/cm  and  Te=40K  in  our  calcula- 


c 

tion.  When  ftfc/277= 0.0025  Hz,  we  find  from  (a)  that  £na(t) 
changes  periodically  with  time  in  a  nonsinusoidal  fashion 
due  to  the  nonlinearity  of  Eq.  (5)  with  respect  to  £na(t). 
There  is  a  tt/2  phase  shift  between  £n2L(t)  and  £b(t)  (thick 
solid  curve  and  right  axis)  because  the  former  is  linearly 
driven  by  d£b(t)/dt  instead  of  £b(t )  as  seen  from  Eq.  (5).  As 
ft  6  increases,  the  phase  shift  disappears.  The  peak  of  £na(t) 
first  increases  with  ft^,  but  then  saturates  to  a  maximum 
value  when  ft^/2ir>0.025Hz,  leading  to  the  saturation  of 
the  tunneling-current  peak.  From  Figs.  1(a)  and  1(b)  we  find 
that  CQW[£na(t)]  and  £m(t)  are  out  of  phase  with  each  other. 
This  can  be  easily  understood  from  the  definition  of  £na(0  in 
Eq.  (4).  The  oscillations  of  CQW[^na(r)]/CQW[0]  with  time 
reflect  the  switching  between  the  charging  and  discharging 
status  of  the  QWs.  There  is  an  asymmetry  in  the  decay  time 
constants  due  to  different  ratios  of  CqW[£na(f)]/CQW[0] 
when  the  QWs  are  in  the  charged  (£na(0<0)  or  discharged 
(£na(0>0)  Status.4 

In  Fig.  2  we  present  the  calculated  total  current 
£/> W  +  £na( t)>Te]  +  €0 erSd £b(t)fdt  [left  axis  in  (a)]  and 
measured  signal  (proportional  to  the  total  current)  [left  axis 
in  (b)]  as  a  function  of  £tbt!2rr  for  different  frequencies 
ftfc/27r.  In  Fig.  2(a),  we  have  picked  7T*=40K  and  £0 

—  1  kV/cm,  and  the  thin  solid,  dashed,  dotted,  dash-dotted, 
and  dash-dot-dotted  curves  correspond  to  ft ^/2  77= 0.02,  0.2, 
0.6,  1,  and  2  Hz,  respectively.  In  Fig.  2(b),  r*=40K  and 
Vo  =  9  V  are  chosen,  and  the  thin  solid,  dashed,  dotted,  dash- 
dotted,  and  dash-dot-dotted  curves  represent  the  results  with 
ft&/27r=50,  20,  10,  5,  and  2  kHz.  £b(t)  and  the  ac  bias 
voltage  are  also  plotted  for  comparison  using  the  thick  solid 
curves  in  (a)  and  (b)  (right  axis).  As  can  be  seen  from  (a),  the 
total  current  is  almost  in-phase  with  £b(t)  when  Slbl2TT 
=  0.02  Hz.  However,  a  tt/2  phase  shift  is  reached  between 
the  total  current  and  £b(t)  when  ft^  becomes  large.  The 
strength  of  the  shifted  peak  increases  with  ft  ^ .  This  behavior 
can  be  explained  as  a  result  of  the  increased  importance  of 
the  dielectric  current  compared  to  the  tunneling  current  as 
ftfc  increases,  as  discussed  in  Sec.  II B.  This  explanation  is 
qualitatively  supported  by  our  experimental  results  shown  in 
(b),  where  the  peak  position  of  the  measured  total  current  is 
gradually  shifted  from  in-phase  to  out-of  phase  compared  to 
the  applied  bias  voltage  and  the  peak  strength  is  enhanced 
with  the  increase  of  ft* . 

Figure  3  displays  the  extracted  tunneling  current  and 
dielectric-current  peaks  [in  (a)]  from  our  calculation  and 
total-current  peaks  [in  (b)]  from  the  measurement  as  a  func¬ 
tion  of  £lb!2ir  on  a  logarithmic  scale.  In  Fig.  3(a)  the  thick 
solid,  dashed,  and  dash-dotted  curves  correspond  to  the 
tunneling-current  peaks  at  7’<?  =  40K  and  £0=  10  kV/cm,  Te 

-  10K  and  £0~  10  kV/cm,  and  Te  =  40  K  and  £0 
=  1  kV/cm,  respectively.  The  thin  solid  and  dash-dotted 
curves  represent  the  dielectric-current  peaks  with  £0 
=  10 kV/cm  and  £q=1  kV/cm.  From  Fig.  3(a)  we  find  that 
the  tunneling-current  peak  is  relatively  large  when  ft^  is 
small,  but  the  dielectric-current  peak  eventually  dominates 
when  ft£  is  large  due  to  the  saturation  of  the  tunneling- 
current  peak  with  ft^ .  The  turning  points  ft  *  are  indicated 
by  the  vertical  arrows  in  (a)  for  each  case,  which  depends  on 
Te  and  £0.  When  £lb  goes  beyond  ft£ ,  the  photodetectors 
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FIG.  2.  Calculated  total  current  (sura  of  the  tunneling  and  dielectric  cur¬ 
rents)  from  sample  1  [left  axis  in  (a)]  and  measured  signal  (proportional  to 
the  total  current)  from  sample  2  [left  axis  in  (b)]  at  7^  =  40  K  as  functions  of 
Q.btl2TT  for  different  frequencies  £lbf27T.  £b{t )  and  ac  bias  voltage  are  also 
plotted  for  comparison  using  thick  solid  curves  [right  axis  in  (a)  and  (b)]. 
We  have  set  in  (a)  £0~ 1  kV/cm,  and  the  thin  solid,  dashed,  dotted,  dash- 
dotted,  and  dash-dot-dotted  curves  correspond  to  0^/277=0.02,  0.2,  0.6,  1, 
and  2  Hz.  The  label  X5  indicates  the  factor  of  amplification  for  the  thin 
solid  curve.  In  (b),  we  have  set  V0—9  V,  and  the  thin  solid,  dashed,  dotted, 
dash-dotted,  and  dash-dot-dotted  curves  correspond  to  0^/277=50,  20,  10, 
5,  and  2  kHz.  Here,  the  curves  are  only  shown  for  one  time  period  for 
clarity. 

suffer  from  a  dynamical  breakdown,  can  be  estimated  by 
l/1ZtCg ,  where  Cg—€Q€rSILt  with  total  length  Lt  of  the 
QWIP  and  lZt  is  the  sequential-tunneling  resistance  depend¬ 
ing  on  Te  and  £0 .  If  Te  is  high,  becomes  large  due  to 
lower  1Zt  (by  comparing  dashed  and  solid  arrows).  For  the 
same  reason,  O*  increases  with  £0  (by  comparing  dash- 
dotted  and  solid  arrows).  This  explanation  is  qualitatively 
confirmed  by  our  experimental  results  shown  in  Fig.  3(b), 
where  the  solid  and  dashed  curves  represent  the  results  at 
Te=40  and  10  K,  and  the  symbols  (□,  A,  ★,  and  •)  corre¬ 
spond  to  the  total-current  peaks  with  V0=  10,  7.5,  5,  and  3  V, 
respectively. 

B.  Emission-current  spike 

For  the  time-dependent  incident  optical  flux  <J)ph(0*  we 
will  use  a  periodic  or  a  step  function  for  an  optical  chopper 
or  an  optical  shutter,  respectively.  4>ph(f)  is  given  by 


Qb/2n  ( kHz ) 

FIG.  3.  Calculated  logarithmic  tunneling-current  (thick  curves)  and 
dielectric-current  (thin  curves)  peaks  from  sample  1  [in  (a)]  and  measured 
logarithmic  total-current  peaks  from  sample  2  [in  (b)]  as  functions  of 
Clbtl2Tr  in  a  logarithmic  scale  for  different  values  of  Te  and  £0  in  (a)  or  Te 
and  V0  in  (b).  In  (a),  the  thick  solid,  dashed,  and  dash-dotted  curves  repre¬ 
sent  the  tunneling-current  peaks  with  7^=40  K  and  £^=10  kV/cm,  Te 
=  10  K  and  £0  =  10 kV/cm,  and  Te= 40 K  and  £0=1  kV/cm,  respectively. 
The  thin  solid  and  dash-dotted  curves  correspond  to  the  dielectric-current 
peaks  with  €0=  10  and  1  kV/cm.  The  solid  and  dashed  curves  in  (b)  repre¬ 
sent  the  results  of  measured  total-current  peaks  at  7^  =  40  and  10  K,  respec¬ 
tively,  with  V0=  10  (□),  7.5  (A),  5  (★),  and  3  V  (•).  The  vertical  arrows  in 
(a)  and  (b)  indicate  the  turning  points  of  f l b  ,  and  the  dynamical  breakdown 
occurs  beyond  those  values. 


f  4>o  +  A<&P(f)  for  an  optical  chopper 
*>h(f)  [Oo+A$0(r)  for  an  optical  shutter  * 

(11) 

where  A<f>  is  the  step  height  and  <f>0  is  the  background  opti¬ 
cal  flux.  Moreover,  6{t)  is  a  step  function  which  steps  up  at 
f=0,  and  P(t)  is  a  unit  periodic  square- wave  function  given 
by 


0  if  mTp^t<(m+V2)Tp 
1  if  (m+l/2)Tp^t<(m  +  l)T' 


In  Eq.  (12),  m  =  0,  ±1,  ±2,...,  Tp  is  the  time  period,  and 
P(t )  steps  up  at  t=(m  +  l/2)Tp  and  steps  down  at  t=(m 
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FIG.  4.  Calculated  £na(r)  [in  (a)]  and  transient  responsivity  Rsp(t)  [in  (b)] 
from  sample  1  at  Te= 40  K,  £b= 2  kV/cm,  and  1 X 1011  cm^s"1  as 
functions  of  time  for  different  step  heights  A<I>=1  (thin  solid  curves),  5  (thin 
dashed  curves),  and  10Xl015cm-2s-1  (thin  dotted  curves)  in  the  presence 
of  an  optical  shutter  opening  at  t  —  0.  The  static  responsivity  is  also  shown 
by  a  thick  dashed  line  in  (b)  for  comparison. 


+  1  )Tp.  In  our  numerical  calculation,  the  excited-state  life¬ 
time  Ti  in  Eq.  (8)  is  taken  to  be  1  ps,  the  electron  temperature 
is  Te  —  40  K,  and  a  dc  electric  field  £b=2  kV/cm  is  applied  to 
the  sample  in  the  z  direction.  In  the  following,  we  still  use 
sample  1  in  our  numerical  calculation.  Some  experimental 
results  from  Ref.  2  are  also  shown  for  a  qualitative  compari¬ 
son  with  our  calculated  results. 

Figure  4  displays  the  solution  of  Eq.  (8)  in  the  presence 
of  an  optical  shutter  which  opens  at  ?= 0,  where  the  nona- 
diabatic  electric  field  £m(t)  [in  (a)]  and  the  transient  respon¬ 
sivity  Rsp(t)  [in  (b)]  are  shown  as  a  function  of  time  for 
different  step  heights  A<t>  =  lX1015  (thin  solid  curves),  5 
X1015  (thin  dashed  curves),  and  10X  10I5cm“2s_l  (thin 
dotted  curves).  Here,  we  have  chosen  £b=  10 kV/cm,  Te 
=  40K,  and  <F0  =  1  X 1011  cm"2s_1  in  our  calculation.  The 
static  responsivity  is  also  shown  by  the  thick  dashed  line  in 
(b)  for  comparison.  From  (a)  we  know  that  £nd(t)>0  after 
the  shutter  is  opened,  indicating  a  discharged  status  for  the 


• 

QWs  during  a  series  of  intermediate  transient  states.  £m(t) 
has  an  initial  rise  just  after  the  shutter  is  opened.  Once 
d<£>ph(t)/dt=0  is  reached,  £na(f)  goes  through  a  decay  down 
to  zero.  Consequently,  the  electrons  in  the  QWs  are  stabi¬ 
lized  by  entering  into  a  final  steady  state  with  a  higher  emis¬ 
sion  current  but  the  same  tunneling  current.  The  initial  rise  of 
£na(0  becomes  steeper  and  steeper  and  the  final  decay  be¬ 
comes  faster  and  faster  as  A3>  is  increased.  From  (b)  we  find 
that  transient  responsivity  Rsp(t)  is  always  smaller  than  its 
static  counterpart  because  of  the  loss  of  electrons  in  the  dis¬ 
charged  QWs.  The  initial  drop  of  Rsp(t)  due  to  the  steep  rise 
of  £na(0  becomes  deeper  and  deeper,  which  is  followed  by  a 
final  decay  to  its  static  value  with  an  enhanced  rate  as  A4> 
increases.  In  the  adiabatic  case,  the  photoemission  current  is 
proportional  to  <bph(0,  which  gives  rise  to  the  static 
responsivity.5,6  The  transient  responsivity  depends  on  the 
time  through  £na(t)  which  depends  on  <l>ph(0  in  the  nona- 
diabatic  case,  implying  a  nonlinear  relation  between  the  pho¬ 
toemission  current  and  d>ph(£). 

In  Fig.  5  we  present  the  calculated  total  ac  cur¬ 
rent  /ff[^+5oa(0,®ph(0]+A/,(0  with  A /,(*) -/,[£* 
+ ^na(0 ’Te\— It\_ £b , Te]  [in  (a)]  and  the  measured  total  ac 
signal  (proportional  to  the  total  ac  current)  [in  (b)]  for  an 
optical  shutter  as  a  function  of  time.  In  (a),  we  set  £b 
=  2  kV/cm,  T,=40K,  and  <I>0=  1 X 1011  cm“2s_1  in  our 
calculation  and  the  solid,  dashed,  dotted,  dash-dotted,  and 
dash-dot-dotted  curves  correspond  to  the  step  heights 
A<&=1,  5,  10,  20,  and  40X1015cm“2s_1,  respectively.  In 
(b),  a  dc  bias  voltage  Vb~  1  V  is  applied  and  — 43  K  and 
4>0=l.lX1013cm“2s“1  are  chosen.  The  11  curves  from  the 
bottom  to  the  top  represent  the  results  with  Ad>=4.5,  6.1, 
7.0,  8.07,  8.5,  12.9,  15.9,  18.9,  21.9,  and  25.9 

X 1013  cm-2  s"1,  respectively.  From  (a),  we  find  that  the  step 
height  of  the  total  ac  current  increases  with  A<&  after  the 
shutter  is  opened.  The  rise  and  decay  time  becomes  shorter 
as  A4>  increases.  Moreover,  an  emission-current  spike  occurs 
when  A$  =  4X  lO^cnT^s"*1.  The  decreasing  rise  and  de¬ 
cay  time  results  from  the  reduced  CQW[£na(0]  due  to  the 
enhancement  of  a  positive  £na(0  with  Ad>  as  shown  in  Fig. 
4(a).  The  emission-current  spike  is  a  result  of  the  competi¬ 
tion  between  the  exponential  enhancement  of  the  escape 
probability  with  respect  to  £na(0  and  the  linear  reduction  of 
electron  density  with  respect  to  £na(0  due  to  the  loss  of 
electrons  in  each  QW,  as  explained  in  Sec.  II C.  This  expla¬ 
nation  is  qualitatively  supported  by  the  experimental  data 
from  Ref.  2  in  Fig.  5(b),  where  both  the  occurrence  of  the 
emission-current  spike  and  the  shortening  of  the  rise  and 
decay  time  are  observed  with  increasing  A<I>. 

Calculated  total  ac  currents  Ie[£b+£ na(0^ph(0] 
+  A It(t)  (thick  solid  curves)  with  an  optical  chopper  are 
shown  in  Fig.  6  as  a  function  of  time  for  different  time  pe¬ 
riods  7^  =  0.1  [in  (a)]  and  0.5  s  [in  (b)].  Here,  we  set  £b 
=  2  kV/cm,  Te=40K,  3>0=1  X  10n  cm‘2s_1,  and  A3> 
=  4.3  X  1016cm_2s_1  in  our  calculation.  The  adiabatic  emis¬ 
sion  currents  Ie[£b  ,Oph(0]  are  also  shown  in  both  Figs.  6(a) 
and  6(b)  by  the  thin  dotted  curves  for  comparison.  From  (a) 
and  (b),  we  find  that  the  rising  time  is  always  relatively 
shorter  than  the  descending  time  within  the  decaying  region 
with  The  rising  and  descending  processes  are 
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FIG.  5.  Calculated  total  ac  current  Ie[€b+£M(t),  $ph(f)3  +  A/,(0  from 
sample  1  [in  (a)]  with  transient  part  of  tunneling  current  A It(t)—It\_£b 
+  £nz{t) ,T g]- 1 :[£b  ,Te]  and  the  measured  total  ac  signal  (proportional  to 
the  total  ac  current)  [in  (b)]  from  Ref.  2  as  functions  of  time  for  an  optical 
shutter  opening  at  /=0.  We  set  in  (a)  £b=2  kV/cm,  Te=40  K,  and  $0=1 
X10ncm“2s-1  in  our  calculation,  and  the  solid,  dashed,  dotted,  dash- 
dotted,  and  dash-dot-dotted  curves  correspond  to  A<E>=1,  5,  10,  20,  and  40 
X 1015  cm-2  s'1,  respectively.  In  (b),  a  dc  bias  voltage  Vb=  1  V  is  applied  to 
the  sample,' and  Te= 43  K  and  <I>0=  1.1X  1013cm~2s~1  are  chosen.  The  11 
curves  from  the  bottom  to  the  top  represent  the  results  with  A<J>=4.5,  6.1, 
7.0,  8.07,  8.5,  12.9,  15.9,  18.9,  21.9,  and  25.9X  1013  cm-2s_1,  respectively. 
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FIG.  6.  Calculated  total  ac  current  /e[£*+£ltt(0,<&ph(0]  +  A/,(0  (thick 
solid  curves)  with  A It(t)=It[£b+£m(t),Te]-It[£b  ,Te]  from  sample  1 
at  Te=40K,  £*,=2  kV/cm,  <l>0=  1 X 1011  cm“2s~\  and  A  <£=4.3 

X  1016cm-2s-1  as  a  function  of  time  for  an  optical  chopper  with  Tp =0.1 
[in  (a)]  and  0.5  s  [in  (b)].  For  comparison,  the  adiabatic  emission  currents 
Ie[£b  ,<£ph(0]  are  also  shown  by  the  thin  dotted  curves  in  both  (a)  and  (b). 
Here,  the  curves  are  only  shown  for  two  time  periods  for  clarity. 


slower  in  (a)  with  respect  to  (b).  A  larger  residual  ac  current 
is  found  in  (a)  at  the  end  of  a  decaying  process.  These  fea¬ 
tures  predicted  by  our  theory  agree  well  with  the  measured 
total  ac  signals  from  Ref.  2  presented  in  Fig.  7  for  Tp=5  [in 
(a)]  and  20  s  [in  (b)],  where  Vb=l  V,  Te= 43  K,  <1>0=1.1 
X1013cm"2s-1,  and  A<I>  =  6. 1 X 1013  ernes'"1  are  chosen 
in  the  experiment. 


C.  Thermal  hysteresis  loop 

For  the  time-dependent  device  temperature,  we  pick  a 
broadened  staircase  form 


re(/)=r0+Ar2 


i  l 

rH — tan' 


r 


where  j  is  the  index  of  the  step,  Ns  is  the  total  number  of 
steps,  ti  is  the  first  stepping  time,  T0  is  the  initial  tempera¬ 
ture,  and  At,  AT,  and  T  are,  respectively,  the  delay  time, 
step  height,  and  step  broadening.  In  our  numerical  calcula¬ 
tion,  the  energy-relaxation  time  rr  in  Eq.  (3)  is  taken  to  be  1 
ps.  The  temperature  evolution  with  time  can  be  realized  by 
turning  on/off  the  heater  at  the  stepping  time  for  stepping 
up/down,  respectively.  In  the  following,  we  choose  sample  1 
for  our  numerical  calculation. 

Figure  8  shows  the  calculated  nonadiabatic  (left  axis  and 
solid  curves)  and  adiabatic  (left  axis  and  dashed  curves)  tun¬ 
neling  currents  with  Sb~5  kV/cm,  Ns—  1,  f1  =  500s  and 
T— 30  s  as  a  function  of  time  for  different  initial  tempera¬ 
tures  T0  [in  (a)]  and  for  positive/negative  AT  [in  (b)].  In  (a), 
AT=5  K  is  chosen,  and  the  thick  and  thin  curves  represent 
the  results  with  T0  =  30  and  10  K.  The  time-dependent  device 
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.  FIG.  7.  Measured  total  ac  signals  (proportional  to  the  total  ac  currents)  from 
Ref.  2  as  functions  of  time  for  an  optical  chopper  with  Tp  =  5  [in  (a)]  and  20 
s  [in  (b)].  In  the  experiment,  a  dc  bias  voltage  Vb=\  V  is  applied  to 
the  sample,  and  7^=43  K,  1.1X  1013  cm"2  s"1,  and  Ad>  =  6.1 
X  1013  cm  2  s  1  are  chosen.  Here,  the  curves  are  only  shown  for  one  time 
period  for  clarity. 


temperature  Te{t )  is  also  displayed  by  the  dash-dotted  curves 
(right  axis).  In  (b),  T0=  10 K  is  taken,  and  the  thick  and  thin 
curves  correspond  to  AT=  -  5  and  5  K.  From  Fig.  8(a),  we 
find  that  the  nonadiabatic  current  (solid  curves)  is  always 
smaller  than  the  adiabatic  one  (dashed  curves)  when  the  tem¬ 
perature  steps  up.  The  amplitude  of  the  reduction  of  It[£b 
+  £na(t)Je(t)]  with  respect  to  It[£b,Te(t)]  increases  with 
decreased  T0,  and  gradually  develops  into  a  dip  when  T0  is 
down  to  10  K.  The  reason  for  this  feature  is  that  charge 
accumulation  (£na(t)<0)  in  the  QWs  blocks  the  sequential 
tunneling  of  electrons  into  the  QWs,  as  explained  in  Sec. 
HA.  From  Fig.  8(b),  we  see  I,[Sb+£J<t),Te{t)'\ 
re(0]  (thin  curves)  when  Te(t)  steps  up,  as  dis¬ 
cussed  in  (a),  but  It[€b+£na(i),Te(t)]>I,[£b,Te(t)]  (thick 
curves)  when  Te(t)  steps  down.  This  implies  a  counterclock¬ 
wise  thermal  hysteresis  loop  in  the  tunneling  current  as  a 
function  of  Te(t),  which  will  be  further  addressed  below. 

Figure  9  presents  the  solutions  £m(t)  of  Eq.  (3)  [left  axis 
and  thick  curves  in  (a)]  and  the  calculated  nonadiabatic  tun¬ 
neling  currents  I,[£b  +  £JJ),Te{t)]  [thick  curves  in  (b)]  as 
functions  of  time  with  £b=\  kV/cm,  Ns- 8,  ^  =  200 s, 
r=30  s,  and  A t=  100  s  for  A T=  —5  K  (solid  curves)  and  5 
K  (dashed  curves).  For  convenience,  the  device  temperature 
profiles  Te(t)  [right  axis  in  (a)]  and  the  calculated  adiabatic 
tunneling  currents  I,[£b  ,Te{t)~\  [in  (b)]  as  functions  of  time 
are  also  displayed  for  A  T=  -  5  K  (thin  solid  curves)  and  5  K 
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FIG.  8.  Calculated  nonadiabatic  tunneling  currents  lHSb+S^{t),Te{t)'\  (left 
axis  and  solid  curves)  and  adiabatic  tunneling  currents  Zr[5b  ,Te(t)]  (left 
axis  and  dashed  curves)  from  sample  1  at  Sb~5  kV/cm  as  functions  of  time 
for  different  initial  temperatures  T0  [in  (a)]  and  positive/negative  step 
heights  A T  [in  (b)],  where  A^=l,  T=30  s,  and  t{  =500  s.  We  have  set  in  (a) 
A T~  5  K,  and  the  thick  and  thin  curves  correspond  to  T0=  30  and  10  K.  For 
comparison,  Te(t )  as  a  function  of  time  is  also  shown  by  the  dash-dotted 
curves  (right  axis).  In  (b),  T0=  10  K  is  chosen,  and  the  thick  and  thin  curves 
represent  the  results  with  A  T——5  and  5  K. 


(thin  dashed  curves).  From  Fig.  9(a),  we  find  that  £na(r) 
(thick  solid  curve)  remains  positive  and  exhibits  an  oscillat¬ 
ing  feature  when  the  device  temperature  is  swept  down  from 
50  to  30  K,  indicating  a  modulated  discharging  status  for  the 
QWs.  On  the  other  hand,  when  the  device  temperature  is 
swept  up  from  30  to  50  K,  £na(Q  (the  thick  dashed  curve)  is 
found  to  be  negative  and  oscillating  with  time,  resulting  from 
a  modulated  charging  status  for  the  QWs.  From  Fig.  9(b),  we 
see  /r[£fc  +  £na(0,rc(r)]  (thick  solid  curve)  shows  a  strong 
oscillation  and  an  enhancement  compared  to  It[£b  ,Te(t)] 
(thin  solid  curve)  when  the  temperature  is  swept  down.  The 
reason  for  this  is  that  the  loss  of  electrons  in  the  QWs  in  a 
modulated  discharging  status  promotes  the  sequential  tunnel¬ 
ing  of  electrons  into  the  QWs  through  the  space-charge-field 
effect.4  Based  on  the  same  concept  of  a  space-charge-field, 
we  can  easily  understand  that  It[£b+ (thick 
dashed  curve)  will  have  a  partial  suppression  compared  to 
h\.£b  >re(f)]  (thin  dashed  curve)  as  the  temperature  is  swept 
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FIG.  9.  Calculated  nonadiabatic  electric  field  £na(  t)  (left  axis  and  thick 
curves)  and  profile  of  Te(t )  (right  axis  and  thin  curves)  in  (a)  from  sample 
1,  as  well  as  the  calculated  tunneling  currents  It[£b+£m(t),Te(t)]  (thick 
curves)  and  It[£b  ,Te(t )]  (thin  curves)  in  (b),  at  £b—  1  kV/cm  as  functions  of 
time  for  positive/negative  step  heights  A  7\  where  Ns— 8,  r1  =  200s,  T=30 
s,  and  A t=  100  s.  In  (a),  AT=  —5  and  5  K  correspond  to  solid  and  dashed 
curves.  We  set  in  (b)  A T~  —5  and  5  K  represented  by  the  solid  and  dashed 
curves,  respectively. 


up  due  to  a  blocking  of  the  sequential  tunneling  of  electrons 
by  the  accumulation  of  electrons  in  the  QWs. 

In  Fig.  10  we  show  the  calculated  nonadiabatic  (thick 
curves)  and  adiabatic  (thin  dotted  curves)  tunneling  currents 
using  Ns=  8,  ti  =  200  s,  T=30  s,  and  A/=  100  s  as  a  function 
of  Te(t),  respectively,  for  different  dc  electric  fields  £b 
=  1  kV/cm  [in  (a)]  and  10  kV/cm  [in  (b)]  with  AT-  -5  K 
(thick  solid  curves)  and  5  K  (thick  dashed  curves).  From  Fig. 
10(a),  we  find  that  the  nonadiabatic  tunneling  current  It[£b 
+  £na(0>^(0]  (thick  dashed  curve)  goes  up  with  increased 
Te{t)  accompanied  by  a  partial  suppression  in  comparison 
with  It[£b  ,Te(t )]  (thin  dotted  curve),  as  explained  in  Fig.  9. 
On  the  other  hand,  /*[£*+ £^(0*^(0]  (thick  solid  curve) 
goes  down  with  decreased  Te(t)  accompanied  by  an  en¬ 
hancement  in  comparison  with  It[£b,Te(t)].  As  a  result,  a 
counterclockwise  loop  of  tunneling  current  due  to  the  ther¬ 
mal  hysteresis  loop  is  formed  as  the  device  temperature  is 
first  swept  up  and  then  swept  down.  The  deviation  of  the 
nonadiabatic  tunneling  current  from  the  adiabatic  one  be¬ 
comes  smaller  as  the  device  temperature  is  increased.  This 
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FIG.  10.  Calculated  tunneling  currents  /,[£*+ £M(f)»re(/)]  (thick  curves) 
and  It[£b  ,Te(t)]  (thin  dotted  curves)  from  sample  1  at  £b=  1  kV/cm  as 
functions  of  devive  temperature  Te(t)  for  different  dc  electric  fields  £b—\ 
and  10  kV/cm.  In  both  (a)  and  (b),  the  results  with  positive  and  negative  step 
heights  A T=  ±5  K  are  represented  by  the  thick  dashed  and  solid  curves, 
where  Ns~  8,  fj  =200  s,  T=30  s,  and  At— 100  s.  Here,  the  solid  (down)  and 
dashed  (up)  arrows  indicate  the  directions  that  the  device  temperature  is 
swept  down  or  up. 

thermal  hysteresis  loop  in  the  tunneling  current  was  reported 
before  by  Choi  et  al?  However,  a  different  explanation, 
which  is  attributed  to  the  effects  of  unintentional  dopants  in 
the  QWIP  barriers,  was  provided  for  this  observed  phenom¬ 
enon  compared  with  the  one  given  here.  The  thermal  hyster¬ 
esis  loop  can  be  effectively  suppressed  by  applying  a  larger 
dc  electric  field,  as  shown  in  Fig.  10(b)  because  the  tunneling 
resistance  will  be  reduced  by  increasing  either  the  tempera¬ 
ture  or  the  bias.3 

IV.  CONCLUSIONS  AND  REMARKS 

In  conclusion,  a  nonadiabatic  sequential-tunneling  model 
has  been  developed  and  used  to  explore  the  common  origin 
of  the  transient  behaviors  of  electrons  in  QWIPs  in  the  pres¬ 
ence  of  time-dependent  external  sources  such  as  device  tem¬ 
perature,  electric  field,  and  incident  optical  flux.  The  nona- 
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diabatic  effect  of  a  time-dependent  source  includes  the 
charge-density  fluctuation  in  the  QWs  induced  by  a  nonadia- 
batic  current  surge.  The  results  from  our  numerical  calcula¬ 
tions  have  been  qualitatively  compared  with  our  experimen¬ 
tal  data  and  those  from  the  experiment  in  Ref.  2  by  Hubbs 
et  al.  The  experimental  observations  have  been  reproduced 
numerically  and  explained  physically. 

For  a  time-dependent  device  temperature,  a  counter¬ 
clockwise  hysteresis  loop  for  the  tunneling  current  as  a  func¬ 
tion  of  swept  temperature  has  been  predicted  and  attributed 
to  the  blockade  or  promotion  of  sequential  tunneling  of  elec¬ 
trons  into  the  QWs  due  to  a  modulated  space-charge  field 
when  the  device  temperature  is  swept  up  and  then  down. 
When  a  time-dependent  electric  field  is  applied  to  the  sys¬ 
tem,  the  peak  of  the  tunneling  current  was  found  to  saturate 
when  the  oscillation  frequency  of  the  ac  electric  field  was 
large  enough.  The  dynamical  breakdown  of  photodetectors 
has  been  predicted  by  including  a  dielectric  current,  which  is 
proportional  to  the  oscillation  frequency  of  an  ac  electric 
field  and  whose  peak  becomes  larger  than  the  value  of  the 
saturated  tunneling-current  peak  in  the  high-frequency  do¬ 
main.  If  an  incident  optical  flux  becomes  time  dependent,  an 
emission-current  spike  was  predicted  as  a  result  of  the  domi¬ 
nant  enhancement  of  the  escape  probability  of  electrons  in 
QWs  over  the  loss  of  electron  density  when  the  applied  dc 
electric  field  is  small. 

The  report  of  a  clockwise  hysteresis  loop  for  the  emis¬ 
sion  current  as  a  function  of  swept  temperature  by  Choi  et  al 
in  Ref.  3  can  also  be  explained  by  our  model.  From  our 
present  study,  we  find  that  the  QWs  remain  either  in  a 
charged  (£na(0<0)  or  a  discharged  (£na(f)>0)  transient 
status  when  the  device  temperature  is  swept  up  or  down, 
respectively  [see  Fig.  9(a)].  From  our  previous  studies,4  we 
know  that  the  nonadiabatic  tunneling  current  It[£b 
+  £na(0>7V(0]  increases  with  £na(f).  However,  the  nonadia¬ 
batic  emission  current  decreases  as 

£na(t)  increases.5  Consequently,  from  the  counterclockwise 
hysteresis  loop  of  the  tunneling  current  we  would  predict  a 
clockwise  loop  for  the  emission  current3  when  the  device 
temperature  is  swept  up  and  then  down. 

Finally,  we  would  like  to  mention  that  we  have  incorpo¬ 
rated  the  charge -density  fluctuation  into  our  model  simply  by 


deriving  a  dynamical  nonadiabatic  electric  field  which  is  spa¬ 
tially  uniform.  This  nonadiabatic  field  results  from  the  space 
charge  in  the  system,  and  should  be  position  dependent  and 
obtained  from  Poisson’s  equation  in  general.  In  this  way,  the 
doping  profile  inside  the  QWIPs,  in  addition  to  electron  den¬ 
sity,  will  have  an  effect  on  the  magnitude  of  the  nonadiabatic 
field.  Also,  it  has  become  well  known  that  the  imbalance 
between  the  injection  current  through  the  emitter  and  the 
emission  current  from  the  QWs  can  cause  a  nonuniform  dis¬ 
tribution  of  the  applied  bias  voltage  across  the  QWIP 
structure.11,12  However,  this  effect  can  be  minimized  when 
the  number  of  QWs  in  the  QWIP  is  large.  The  present  fea¬ 
tures  predicted  by  our  nonadiabatic  model  in  this  article  re¬ 
main  valid  after  including  the  nonuniform  distribution  of  the 
bias  voltage  inside  the  sample,  although  some  quantitative 
difference  is  expected. 
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